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10.  aBOTJIACT  fCaiMnu.  it  nvm—  mik*  M MMaMir  mk  IOmmIIV  Or  ttoaO  aimkH 

One  of  the  complex  and  unsolved  problems  in  the  sonar  field  is  the 
modeling  of  a mathematical  and  computerized  solution  for  the  maximization 
of  signal-to-background  ratio  (SBR) . In  dealing  with  the  maximization  of  SBR, 
it  has  been  found  sufficient  to  attack  the  problem  of  maximizing  signal-to- 
reverberation  ratio  (SRR) . In  this  report,  a numerical  mathematical  model 
is  introduced  for  the  complete  solution  of  the  SRR-maximization  problem.  The 
techniques  employed  to  develop  the  entire  optimization  process  are  ^“l|)^(over) 
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0.  ABSTRACT  (Continued) 


justified  by  mathematical  theory.  A mathematical  expression  for  SRR  is 
obtained,  the  complete  SRR-optimization  problem  is  addressed,  and  the 
development  of  a systematic  computerized  numerical  procedure  to  optimize 
SRR  is  outlined.  Also,  the  probleau  encountered  in  the  mathematical  and 
computer  modeling  are  discussed.  The  computerized  model  has  been  applied 
to  several  sonar  problems.  Numerical  results  have  demonstrated  the  validity 
of  the  mathematical  model  and  the  effectiveness  of  the  computer  BK>del. 
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MODELING  A MATHEMATICAL  SOLUTION  FOR  THE 
OPTIMIZATION  OF  SIGNAL-TO-REVERBERATION  RATIO 


INTRODUCTION 


A sonar  systea  uses  acoustic  energy  for  the  purposes  of  detection, 
distance  estiaation,  and  coanuni cations.  A passive  sonar  receives  acous> 
tic  energy  generated  by  a target,  whereas  an  active  sonar  generates 
acoustic  energy  that  is  reflected  by  a target  and  detected  by  the  sonar 
1 receiver, 

i 

1 There  is  a desired  and  an  undesired  acoustic  energy  associated  with 

[ the  natural  acoustical  environment.  The  desired  portion  is  the  signal 

and  the  undesired  portion  is  the  background.  The  background  is  either 
ambient  noise  or  reverberation  noise.  The  main  objective  of  research  is 
to  increase  sonar  performance,  which  means  increasing  the  overall 
j response  of  the  sonar  system  to  the  signal  and  decreasing  the  overall 

response  of  the  sonar  system  to  the  background,  i.e.,  increasing  the 
I signal-to-background  ratio  (SBR). 

I . 

i One  of  the  complex  and  unsolved  sonar  problems  is  the  modeling  of 

a computerized  mathematical  solution  for  the  optimization  of  the  SBR. 
Since  optimization  of  signal-to-noise  ratio  (SNR)  is  a special  case  of 
the  optimization  of  signal-to-reverberation  ratio  (SRR) , it  is  sufficient 
to  consider  the  general  SRR  problem  in  dealing  with  the  optimization  of 
SBR.  Much  of  the  existing  theory  and  techniques  that  have  been  developed 
have  been  directed  at  maximizing  SNR,  leaving  unfinished  a comprehensive 
treatment  of  the  general  case.  In  general,  reverberation  produced  by 
the  scatterers  is  not  identical;  therefore,  the  distribution  is  not 
uniform.  This  presents  a formulation  difficulty  so  that  an  exact  mathe- 
\ matical  expression  for  the  SRR  is  not  easy  to  obtain.  However,  under 

I some  realistic  environmental  conditions,  reasonable  assumptions  can  be 

made  in  order  to  make  the  calculation  of  SRR  by  a mathematical  expres- 
I sion  possible.  From  this  mathematical  expression,  a complete  definition 

of  SRR  will  be  established  in  order  to  classify  SBR  in  different  cate- 
gories. The  expression,  a two-dimensional  integral,  is  difficult  to 
evaluate  in  general.  In  many  applications  the  solutions  are  special 
I functions,  such  as  Bessel  functions.  In  others,  a closed-form  mathema- 

tical solution  is  extremely  difficult  to  obtain.  Reverberation  power 
received  from  the  array  is  positive;  this  physical  property  allows  the 
application  of  the  quadratic- form  theory  to  optimize  SRR.  Note  that 
other  techniques  used  to  maximize  SNR,  such  as  Lagrange  Multipliers,^ 
matching  filter  technique, 2 statistical  techniques, 3 and  others, are 
all  applicable  to  maximize  SRR.  However,  this  model  uses  the  elegant 

I 
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quadratic- form  theory  because  it  yields  an  exact  fomulation  for  the 
systematic  determination  of  signal  and  reverberation  matrices  that  is 
very  suitable  for  any  computer.  Under  general  conditions,  when  the 
mathematical  expression  cannot  be  obtained,  background  measurements  can 
be  used  together  with  the  quadratic- form  theory  to  optimize  SRR. 

A conqjlete  computerized  procedure  that  implements  the  theory  will 
be  outlined  in  this  report.  A detailed  description  of  the  computer 
model  will  be  reported  separately.  A natural  question  arises:  Can  this 
combined  procedure  be  improved?  A discussion  of  the  problems  encoun- 
tered and  some  practical  solutions  may  serve  as  an  answer  to  the  above 
question.  A sizeable  cylindrical  array  is  presented  as  an  exaiq>le  to 
illustrate  the  problems  encountered  and  to  show  how  this  problem  was 
solved  practically.  Experimental  computer  programs,  designed  to  imple- 
ment this  procedure,  are  written  in  ANSI  FORTRAN  language,  but  will  not 
be  included  in  this  report.  However,  the  finalized  package  of  programs 
will  be  included  in  the  cooqjuter  model  report. 


A MATHEMATICAL  MODEL  TO  MAXIMIZE  SRR 


In  this  report,  dj , Rj,  and  R are  functions  of  0,^;  U,  V are  func 
tions  of  e,(>,  and  aj.  "^For  economy  in  wording,  we  have  dropped  the 
symbols  6,^,  and  aj  for  the  above  functions. 

BEAM  OUTPUT  FUNCTION  AND  PHYSICAL  CONSIDERATIONS 

Figure  1 shows  a configuration  in  three-dimensional  space  with 
spherical  angles  0,4),  which  have  the  following  definitions: 


z 


Figure  1.  Three-Dimensional  Space  With  Spherical  Angles  Q 
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let  6(0  £ 6 ^ 2ii)  define  the  aziauthal  deviation  in  the  horizontal 
plane; 

let  ^ %")  define  the  vertical  deviation  fro«  a horizontal 

reference; 

let  Rj  (0,(»  (or  Rj)  be  the  anplitude  response  of  the  j-th  element 
of  a sonar  array  consisting  of  n elements; 

let  dj  (6,6)  (oi*  dj)  be  the  distance  from  the  j-th  element  to  the 
reference  plane  perpendicular  to  the  direction  (6,6)  of  the  incoming 
acoustic  wave;  and 

let  X be  the  wavelength. 

Distances  are  obtained  by  means  of  direction  cosines  that  define 
the  direction  ot  propagation  perpendicular  to  the  incoming  wavefront. 

The  beam  output  function  V(6,6,arj)  (or  V)  is  defined  by 

V = y a.R.e  y 

^ J J (1) 

where  Oj  are  scalars  to  be  determined. 

Associated  with  each  sonar  system  there  is  a beam  pattern  that  is 
calculated  by  means  of  the  beam  output  function.  An  inportant  sonar 
parameter  for  sound  reception  used  to  measure  the  amount  of  response 
through  the  beam  pattern  and  to  discriminate  against  noise  in  favor  of 
signal  is  conventionally  termed  directivity  index  (DI).  DI  refers  to  a 
plane-wave  signal  in  an  isotropic  noise  field. These  physical  condi- 
tions permit  DI  to  be  expressed  mathematically  by  the  two-dimensional 
integral  given  by  equation  (6) . When  signal  and  noise  fields  have  other 
directionalities,  the  term  "array  gain"(AG)  is  used.  The  mathematical 
expression  of  AG  is  given  by  equation  (5).  The  definition  of  AG  is, 
then,  generalized  to  obtain  a mathematical  expression  for  reverberation 
gain  (RG)  given  by  equation  (3). 

MATHEMATICAL  EXPRESSIONS 

In  the  region  0 £ 0 £ 2Tr,  -ijir  ^ a £ 6 b £ a generalized  direc- 
tive gain,(D)  can  be  defined  as  follows: 


2iT(sin  b - sin  a) 


a 


Ml 

|Uj2 


cos  6d6d6 
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where 

R ■ R(6,^)  for  an  active  sonar  systea  defines  a reverberation 
function  that  is  used  to  characterize  the  relative  reverberation  strength 
of  backscattered  energy  fron  boundary  or  voluae  scatterers  over  the 
angular  interval  of  interest,  and  for  a passive  sonar  systea  defines  a 
noise  function; 

U > U(6,^)  defines  a transmitting'beaa  output  function; 

Uq  » U(6o,(^q)  defines  a transmitting  steering-direction  vector; 

V a V(8,^,aj)  defines  a receiving-beam  output  function;  and 

Vo  = V(eo,(>o,or.)  defines  a receiving  steering-direction  vector. 

J • 

In  equation  (2),  let  D be  D;  when  R > 1:  let  D be  Dj  when  U ■ 1, 
a = -Jjii,  b = let  D be  D3  when  R>1,  U>1,  a«  -Jjir,  b • >jir.  The 

terms  reverberation  gain  RG,  reverberation  index  RI,  array  gain  AG, 
and  directivity  index  DI  then  have  the  following  mathematical  defini- 


tions : 

RG  (reverberation  gain)  » 10  log^gD,  (3) 
RI  (reverberation  index)  * 10  log^oDii  (4) 
AG  (array  gain)  = 10  logioD2,  and  (S) 
DI  (directivity  index)  = 10  logxoD3.  (6) 


It  is  seen  that  RG  in  equation  (3)  is  the  most  general  expression, 
and  that  RI,  AG,  and  DI  are  special  cases  of  RG.  This  shows  that  in 
dealing  with  the  SBR-maximization  problem  it  is  sufficient  to  deal  with 
the  maximization  of  equation  (3),  which  is  equivalent  to  the  maximiza- 
tion of  D in  equation  (2).  Equations  (3)  and  (4)  are  valid  for  active 
sonars  with  reverberation  background,  and  equations  (5)  and  (6)  are 
valid  for  passive  sonars  with  noise  background. 

The  two-dimensional  integral,  occurring  in  equation  (2),  is  usually 
not  simple  to  evaluate.  In  some  applications,  the  exact  solutions  exist 
in  closed  forms;  in  this  case,  the  expression  is  usually  a combination 
of  Bessel  functions. Vfhen  the  closed-form  solutions  cannot  be  obtained, 
numerical  integration  techniques  are  used. 
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THE  SRR-MAXIMIZATION  PROBLEM  AND  A FEASIBLE  SOLUTION 

Under  all  conditions,  it  is  desirable  to  increase  the  signal'to- 
background  ratio  (SBR) . For  the  situations  where  aobient  noise  is  the 
dominant  interference,  the  maximization  of  SNR  is  required.  A closely 
related  and  more  general  problem  arises  when  interference  is  dominated 
by  reverberation;  in  this  situation  the  maximization  of  SRR  is  required. 
In  dealing  with  the  maximization  of  SBR,  it  is  sufficient  to  deal  with 
the  maximization  of  SRR. 

The  SRR-maximization  problem  can  be  described  by  one  question:  Can 
a set  of  scalar  coefficients  aj  be  determined  such  that  the  generalized 
directive  gain  D of  equation  (2)  is  at  its  maximum?  The  answer  is  yes, 
if  D possesses  the  necessary  properties.  The  reverberation  power 
received  by  the  array  is  always  positive;  this  implies  that  the  denomi- 
nator of  D is  positive.  This  property  allows  D to  be  expressed  as  an 
SRR  and  suggests  that  D can  be  expressed  simultaneously  as  a ratio  of 
two  quadratic  forms  with  a positive-definite  denominator;  therefore, 
the  elegant  quadratic- form  theory  can  be  applied  to  maximize  SRR.  The 
quadratic- form  theory  not  only  gives  the  exact  mathematical  formulation 
and  the  optimum  solution  to  the  SRR-maximization  problem,  but  also  gives 
a systematic  optimization  procedure  that  is  very  suitable  for  inple- 
mentation  on  large-scale  computers.  The  approach  is  to  express 
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for  j > A simple  verification  shows  that 

If  - A . (10) 

Equation  (10)  gives  an  alternate  determination  of  aj^.  Similarly, 
setting  equal  to  the  denominator  of  D,  the  n^  elements  of  B ■ (bjg) 

are  expressed  by 


jk  2ii(sln  b - sin  a) 


Cl  ■ r 


lul^  if(cl,-d) 

' * r e cos  ♦ 


d^de.  (11) 


where  ar*Ba  is  positive  because  of  the  physical  property;  this  implies 
B is  positive-definite.  Then  D can  be  maximized  according  to  the 
following  matrix  theory. 

If  A(or,a)  and  B(a,a)  are  2 quadratic  forms  in  the  variable  aj  and 
Bis  positive-definite,  then  the  characteristic  equation  of  A(a,a)  - 
X B(ar><ir)  is 

p(X)  * det  ( A - XB)  - 0 . (12) 

In  view  of  the  matrix  theory  for  eigenvalue  problems, p(X)  has 
real  roots  Xj  such  that  Xj  ^ Xj+j  for  j » l,2,...,n-l  and  Xj,  X^  repre- 
sent the  bounds  of  equation  (77;  i.e., 

. (15) 

The  right  equality  of  equation  (13)  is  attained  if 

Aa=  \Ba  . (14) 


Then,  it  is  seen  that 


X E*B-1E  >0  . (15) 

n a*Ba 

and  the  eigenvector  «opt  is  the  optimum  vector  that  can  be  calculated 
by 

«opt  = ® (16) 

Equations  (15)  and  (16)  constitute  the  complete  solution  to  the  SRR- 
maximization  problem. 
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Then,  optimum  RG  can  be  calculated  from 

%,  ■ ><>  ‘“Sin  ( ‘‘-opt'  • <*’> 

The  details  of  this  complete  mathematical  development,  with  proofs, 
can  be  found  in  existing  material. 


EFFICIENT  IMPLEMENTATION  OF  THE  SRR  MATHEMATICAL  MODEL 


PROCEDURE  IN  FLOW  CHART 

A large  computer  is  necessary  to  carry  out  numerical  implementation 
of  such  a complex  problem.  The  complete  computational  procedure  is  shown  as 
a flowchart  in  Hgure  2.  Computational  difficulties  have  been  minimized 
in  this  procedure  to  reach  a feasible  solution. 

A set  of  experimental  computer  programs  has  been  developed  in 
FORTRAN  language  at  the  Naval  Underwater  Systems  Center  (NUSC)  in  New 
New  London  on  a Uni vac  1108  to  implement  the  above  procedure.  The 
computer  programs  are  not  yet  fully  automated.  Documentation  of  the 
detailed  computerized  model  will  be  made  available  at  a later  date. 

IMPORTANT  CONSIDERATIONS 

In  order  to  construct  2in  efficient  computer  model,  certain  impor- 
tant considerations  cannot  be  neglected.  These  considerations  may  be 
applicable  for  developing  all  types  of  computer  models.  Several  of 
these  considerations  are  incorporated  into  the  development  of  the  SRR 
computer  model;  in  fact,  these  become  desirable  features  of  this  model. 

A brief  discussion  of  these  features  follows. 

Transportabi lity 

An  effective  computer  model  is  expected  to  be  used  often  and  widely. 
It  is  desirable  to  develop  computer  programs  not  for  one  particular 
machine  but  for  many  different  machines  that  speak  the  same  dialect, 
but  perfectly  port^le  programs  are  not  yet  in  existence.  Nevertheless, 
program  writers  should  attempt  to  develop  progreuns  in  such  a way  that 
the  least  amount  of  effort  is  required  to  convert  programs  from  one 
machine  to  the  other.  The  SRR  programs  are  designed  in  ANSI  FORTRAN, 
which  possesses  this  distinguishing  quality,  transportability. 

Flexibility 

Particular  portions  of  the  SRR  programs  are  treated  by  special  cal- 
culations. These  calculations  are  performed  either  by  means  of  existing 

i 
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library  packaged  programs  or  by  numerical  techniques  developed  for  this 
purpose.  Generally  speaking,  these  programs  may  not  be  the  most  effic- 
ient ones.  The  search  for  more  efficient  algorithms  always  goes  on. 

To  give  certain  flexibility,  these  SRR  techniques  are  designed  modularly 
in  subroutine  formats  so  that  they  can  be  distinguished  and  replaced 
more  easily. 

Accuracy  and  Speed 

In  carrying  out  numerical  integrations,  the  step  size  plays  an 
important  role  that  determines  the  accuracy  and  the  speed;  one  may  be 
sacrificed  to  gain  the  other.  However,  accuracy  is  usually  more  inpor- 
tant  than  speed.  To  meet  less  stringent  accuracy  requirements,  speed 
can  be  improved  by  taking  larger  step  sizes  in  the  numerical  procedures. 
High  speed. should  not  be  acconqilished  by  sacrificing  accuracy  beyond 
certain  reasonable  limits.  One  practical  approach  to  improve  the  speed 
and  still  maintain  good  accuracy  is'  to  specialize  the  applications.  The 
block  Toeplitzl^  reverberation  matrix  is  an  example  of  a special  applica- 
tion for  equally  spaced  arrays. 

Reliability 

It  frequently  happens  that  a user  fails  to  supply  a piece  of  input; 
consequently,  the  program  may  fail  or  give  poor  results.  To  avoid  this 
danger,  all  possible  inputs  are  defaulted  to  reasonable  values  based 
on  physical  knowledge  to  maintain  certain  reliability  in  SRR  programs. 


PROBLEMS  IN  MODELING  AND  SOME  PRACTICAL  SOLUTIONS 


PROBLEMS  ENCOUNTERED 

It  is  natural  that  in  the  process  of  modeling  a solution  for  any 
real-world  problem,  other  problems  may  arise.  Although  the  problems 
encountered  here  are  not  unique,  there  are  no  simple  and  easy  solutions 
to  them.  Some  of  the  problems  encountered  in  modeling  a computerized 
mathematical  solution  to  the  SRR-maximization  problem  are  outlined  below 
in  the  order  of  consideration  in  a question-and- answer  format,  along 
with  some  practical  solutions. 

Q;  The  underwater  acoustic  problem  of  obtaining  a mathematical 
expression  for  SRR  involves  three  types  of  conditions;  environmental, 
tactical,  and  equipment.  Each  condition  involves  a number  of  variables. 
It  is  almost  impossible  to  come  iq>  with  a mathematical  expression  that 
is  general  enough  to  cover  all  these  variables.  How  should  the  general 
problem  be  approached  to  obtain  a viable  solution? 
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A:  A practical  approach  is  to  siaplify  the  problem  in  order  to 
obtain  a reasonable  utheBatical  expression  for  SRR.  If  a sinplified 
|t  version  cannot  be  solved,  the  solution  to  the  general  complicated  problem 

can  never  be  expected.  Levels  of  higher  complexity  may  be  added  after 
the  less  conplex  forms  have  been  accepted  as  valid. 

Q:  How  can  this  problem  be  simplified? 

A;  The  problem  may  be  simplified  by  making  some  reasonable,  realis- 
tic assumptions  so  that  the  whole  physical  problem  can  be  expressed 
mathematically.  A starting  point  for  an  active-passive  sonar  system  of 
n elements  in  a general  environment  would  be  to  consider  the  passive  case 
first,  assuming  omnidirectional  array  elements  and  a plane-wave  signal  in 
an  isotropic  noise  field;  then  the  directivity  index  (DI)  can  be  expres- 
sed mathematically  by  a two-dimensional  integral, given  by  equation 
(6).  From  this,  we  begin  by  generalizing  the  definition  of  directive 
gain  in  order  to  establish  the  counterpart  two-dimensional  integral  of 
DI  for  reverberation  gain,  RG  (given  by  equation  (3)),  for  active  sonar 
systems.  Then  we  demonstrate  that  the  mathematical  expression  of  equa- 
tion (2)  for  RG  was  established  with  sufficient  generality  so  that  RI, 

AG,  and  DI  are  automatically  accommodated. 

Q;  Why  make  these  simple  assunptions? 

A:  The  fact  is  that  if  fewer  simplifications  are  made,  a more  com- 
plicated mathematical  expression  may  result.  By  the  same  token,  if  more 
simplifications  are  made,  a simpler  mathematical  expression  may  result. 
Simple  mathematical  expressions  are  easier  to  handle.  As  long  as  the 
simplified  mathematical  expression  is  physically  meaningful,  this  may  be 
a good  way  to  start  the  solution  of  a complex  problem. 

Q;  Is  this  two-dimensional  expression  mathematically  and  physically 
valid? 

A;  The  existence  of  the  two-dimensional  integral  comes  from  the 
fact  that  the  integrand  is  a bounded,  continuous  function  in  both  6 and 
The  physical  validity  is  confirmed  by  the  actual  experiments  and 
measurements  undertaken  by  sonar  scientists. 

■ 21  problems  in  calculating  the  integral? 

A;  Problems  are  easier  if  the  two-dimensional  integral  can  be 
solved  by  closed-form  solutions,  for  example,  as  in  the  case  of  an  equally 
spaced  line  array  in  an  isotropic  noise  field,  where  the  solution  is  in 

the  form  When  closed-form  solutions  cannot  be  obtained,  numerical 

integration  can  be  used  to  evaluate  the  two-dimensional  integral;  a 
disadvantage  is  that  the  numerical  integration  may  be  very  time-consuming 
if  extreme  accuracy  is  demanded  or  a large  array  is  involved. 
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Q;  Is  the  technique  to  maximize  SRR  in  the  mathematical  model 
theoretically  justified? 

A;  Since  the  physical  property  permits  SRR  to  be  expressible  as  a 
ratio  of  two  quadratic  forms,  elegant  quadratic-form  theory  to  maximize 
SRR  is  applicable. 

Q:  What  are  the  mathematical  problems  involved  in  obtaining  the 
optimum  RG? 

A;  Formulating  a complete  mathematical  procedure  to  obtain  optimum 
RG  requires  a great  deal  of  matrix  algebra  and  a thorough  analysis  of 
this  generalized  eigenvalue  problem.  Following  the  above,  equation  (8) 
has  to  be  worked  out  for  the  determination  of  the  signal  matrix  elements, 
and  equation  CH)  has  to  be  worked  out  for  the  determination  of  rever- 
beration matrix  elements.  Finally,  a solution  to  this  generalized 
eigenvalue  problem  has  to  be  worked  out  to  obtain  equation  (17)  for  the 
calculation  of  optimum  RG. 

Q:  Why  is  a the  optimum  vector? 

A;  Maximizing  SRR  requires  solving  a generalized  eigenvalue  problem 
of  the  form  ( A-  XB)a=  0.  In  the  SRR  model,  because  only  one  steering 
direction  at  a time  is  considered,  the  rank  of  A » 1.  Therefore,  there 
exists  only  one  nonzero,  real  X whose  associated  eigenvector  is  a. 

Q;  What  problems  are  encountered  in  obtaining  by  B“^E  ? 

A;  Since  B is  positive-definite,  B is  invertible.  Solving 
( A-  xB)a  - 0 is  the  same  as  solving  B(B~^A  - XI)  a > 0.  As  long  as 
B is  not  ill-conditioned,  inversion  programs  are  readily  available. 
Because  of  small  array  element  spacings,  ill-conditioned  B may  result, 
in  which  case  an  accuracy  problem  is  anticipated  in  inverting  B . 

Q;  Can  RG  always  be  improved? 

A:  The  improvement  of  RG  depends  mainly  on  the  element  spacings. 
Some  arrays,  e.g.,  an  equally  spaced  line  array  with  s^X  spacing,  are 
already  in  optimum  state,  in  which  case  no  improvement  on  RG  can  be 
expected. 

Q;  The  SRR  mathematical  model  presently  in  practice  is  developed 
to  solve  a restricted  problem.  How  can  the  solution  be  generalized  to 
solve  the  most  general  problem? 

A;  Since  the  reverberation  power  received  by  the  array  is  always 
positive,  the  quadratic- form  theory  is  still  applicable,  lliis  mathema- 
tical model  is  still  valid  for  the  reason  that  a solution  to  the  general 
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problem  can  be  obtained  through  the  use  of  measurements  of  the  noise  or 
reverberation  matrix  and  the  application  of  quadratic* form  theory. 

EXISTING  PROBLEMS 

( 

In  addition  to  the  above  problems,  other  problems  exist  that  will 
be  addressed,  again,  in  a question-and-answer  format.  j 

Q;  How  can  computations  be  done  faster  while  maintaining  reasonable 
accuracy? 

A;  A number  of  problems  may  be  similar,  such  that  they  may  have 
some  complicated  calculations  in  common.  These  complicated  calculations 
need  to  be  done  only  once  and  can  be  saved  for  repeated  uses.  For 
example,  since  the  noise  matrix  B is  not  steering  dependent,  B and/or 
B”^  could  be  saved  for  an  array  being  studied  at  different  steering 
angles. 

Q:  Good  results  depend  on  accurate  inputs.  How  can  the  accuracy  of 
the  inputs  be  ensured? 

A;  This  model  has  no  direct  control  of  the  inputs.  Nevertheless, 
some  inputs  are  defaulted  meaningfully.  The  accuracy  of  the  inputs  is 
the  responsibility  of  the  user. 

Q:  How  can  a large  array  be  handled  when  the  reveiberation  matrix 
Cor  noise  matrix)  exceeds  the  computer  capacity? 

^ In  a number  of  applications  the  Toeplitz  technique  can  be  used, 
requiring  only  that  subsets  of  the  entire  matrix  be  stored  because  of 
repetition  of  the  entries.  When  the  reverberation  matrix  has  Toeplitz 
structure  but  still  exceeds  memory,  auxiliary  storages  are  needed  and 
the  model  has  to  be  modified  to  handle  additional  I/O. 

Q:  What  is  a good  method  to  handle  the  inversion  of  an  ill-condi- 
tioned matrix? 

A;  There  is  no  good  way  to  handle  the  inversion  of  an  ill-condi- 
tioned matrix.  Since,  when  maximizing  SRR,  the  problem  is  to  find  the 
eigenvector  or  of  the  problem  (A  - XB)<irs  0,  the  inversion  of  B can  be 
avoided  if  an  iterative  technique  is  used. 

The  search  for  more  effective  answers  to  the  questions  presented 
still  continues.  Until  effective  solutions  are  found,  the  problems  of 
accurate  inputs,  exceeding  memory  capacity,  speed  and  accuracy,  and 
solving  a very  large  eigensystem  will  remain. 
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AN  ILLUSTRATIVE  EXAMPLE 


The  SRR  model  has  been  applied  to  several  sonar  systems.  In  this 
section  an  example  using  a cylindrical  array  of  192  elements  has  been 
selected  for  presentation.  Computations  were  carried  out  in  single-preci- 
sion arithmetic  on  the  NUSC,  New  London  Laboratory,  Univac  1108  conqputer 
with  an  allowable  memory  capacity  of  65K.  Compilation  was  carried  out 
by  means  of  the  EXEC  8 operating  system. 

This  array  is  selected  for  presentation  for  the  following  multiple 
purposes : 

1.  To  show  the  validity  of  the  SRR  model,  both  in  mathematics  and 
in  computer  mechanism; 

2.  To  exhibit  and  compare  the  nonoptimum  and  optimum  beam  patterns, 
as  well  as  nonoptimum  and  optimum  SNR; 

3.  To  show  how  to  overcome  the  difficulties  and  solve  the  problem 
efficiently,  since  this  sizeable  array  encounters  almost  all  the  problems 
mentioned  in  the  previous  section; 

4.  To  show  that  a general  array  of  192  elements  may  not  be  solv- 
able due  to  the  limit  of  memory  capacity  and  existing  support  softwares 
or  the  impractical  computation  speed.  However,  this  case  has  certain 
properties  that  enable  it  to  be  specialized  to  obtain  a practical 
solution  to  the  problem.  The  solution  speed  will  be  accelerated  and  the 
accuracy  can  still  be  maintained  by  this  specialization;  and 

5.  To  show  how  the  Toeplitz  matrix  can  be  used  to  advantage  in 
this  application. 


ARRAY  STRUCTURE  AND  INPUT  CONDITIONS 

To  make  calculations  on  this  cylindrical  array  of  192  elements, 
some  specific  input  conditions  were  assumed.  The  array  consists  of  a 
total  of  8 rows;  each  row  has  24  elements.  The  radius  of  the  cylinder 
is  8 feet.  The  horizontal  angular  separation  is  5°  between  two  adjacent 
elements  and  the  vertical  spacing  between  elements  is  0.72  foot.  Each 
element  is  assumed  omnidirectional  in  an  isotropic  noise  field.  A plane 
wavefront  is  considered  and  the  array  is  steered  in  the  middle  with  a 
depression  angle  of  45°.  A frequency  of  3500  Hz  is  used. 
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PROBLEMS  INVOLVED 

In  addition  to  minor  problems,  there  are  two  major  difficulties: 
the  time  required  for  coiiq>utation  and  the  limitation  on  conqputer  memory. 

In  general,  for  an  arbitrary  array  of  192  elements,  generation  of 

noise  matrix  elements  will  require  Z ^ ■ 18336  determinations. 

Each  determination  is  the  performance  or  a two-dimensional  numerical 
integration  in  the  region  0 ^ 9 ^ 2t,  -*s!t  £ ♦ ^ ^sit.  For  acceptable 
accuracy,  with  a fixed  mesh  size”Csay  3*  in  each  direction),  7200  inte- 
grand evaluations  will  be  required.  If  each  determination  takes  1 sec- 
ond, total  determination  time  will  be  about  5 hours  and  6 minutes.  The 
additional  calculation  of  nonoptimum  and  optimum  beam  patterns,  SNR,  the 
matrix  manipulation  and  inversion,  and  other  factors,  will  push  the 
total  computation  time  up  to  about  6 hours  to  reach  a solution.  This 
lengthy  computational  time  is  a major  difficulty. 

The  noise  matrix  will  generally  result  in  Hermitian  form.  Handling 
a 192  X 192  matrix  of  coiiq>lex  elements  will  exceed  our  present  computer 
memory  capacity,  which  is  another  major  difficulty. 

A PRACTICAL  SOLUTION 

An  attempt  was  made  to  overcome  the  above  difficulties  by  analyzing 
the  properties  of  this  array.  A specialized  approach  to  the  problem  may  be 
used  in  order  to  obtain  a practical  solution  by  examining  these  properties. 

In  general,  equation  (11)  is  used  to  determine  the  elements  bj](  of 
the  noise  matrix.  A thorough  analysis  reveals  that  for  omnidirectional 
elements  in  an  isotropic  noise  field,  the  denominator  of  equation  (11)  is 
a two-dimensional  integral  that  takes  the  form 


■ 


cos  ^d(^de  . 


(18) 


Since  dj  is  calculated  by  direction  cosines,  the  above  integral  is 
evaluated  in  the  form 


niT  . 2ir 


( Xj^- Xj ) cos  ((is  in  0+ (yj^ -y  j ) cos  (|icos  0- ( Zj^  - z j ) s in  (ji  ] 


cos  ^d^d9. 


(19) 
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A closed-form  solution  for  the  above  integral  has  already  been 
worked  out  and  can  be  found  in  existing  literature. The  solution  is 
given  by 


(2iT)3/2jj^(i2Ls)/(|2-S)’*  , 


C20) 


where 


S = [(x^  - Xj)2  + (yj^  - y.)2  + . z.)^]’*  . 

Since  Jjj(kt)  » (1)**  sin(kt)/(kt)\  equation  (19)  can  be  siii5)lified 
to  the  expression 


The  exact  solution  also  reveals  that  bjj^  are  real;  therefore,  the 
program  can  be  specialized  to  handle  a 192  x 192  symmetric  matrix  within 
our  present  memory  capability.  These  simplifications  also  can  reduce 
the  6 hour  computation  time  to  about  20  minutes. 

NUMERICAL  RESULTS 

Following  the  procedure  in  the  flowchart,  we  first  calculate  the 
nonoptimum  beam  patterns  and  nonoptimum  DI.  The  nonoptimum  horizontal 
beam  pattern  is  given  in  figure  3 as  a solid  line.  The  nonoptimum  DI  is 
foimd  to  be  22.35  dB.  We  then  proceed  to  calculate  the  optimum  DI, 
which  requires  the  determination  of  bjj^,  the  matrix  inversion,  and  the 
calculation  of  optimum  shading  coefficients.  The  optimum  horizontal 
beam  pattern  is  also  shown  in  figure  3,  but  as  a dashed  line.  The 
optimum  shading  coefficients  are  given  in  tables  1 and  2.  The  optimum 
DI  is  found  to  be  22.87  dB,  a gain  of  0.52  dB. 

The  locations  of  the  elements  in  Cartesian  coordinates  (X,Y,Z)  are 
calculated  using  cylindrical  coordinates;  i.e.,  X = r cos  0,  Y = r sin  0, 

Z = Z for  32.5®  £ 6 £ 147.5®  at  a 0-increment  of  5®.  Let  W represent 
either  X,  Y,  or  Z;  and  let  the  indexes  j,  k,  m be  such  that  1 £ j,k  £ 24, 
m = 0,1,..., 7.  One  can  verify  that  Wj+24in  ’ Wk+24m  = Wj  - Wk.  This 
gives  rise  to  a block  Toeplitz  noise  matrix;  each  block  is  a 24  x 24 
matrix.  The  total  determination  of  the  noise  matrix  elements  is  reduced 
to  4584  elements  only,  a saving  of  13,752  locations.  At  the  present  time, 
a block  Toeplitz  inversion  subroutine  is  not  available.  To  handle 
the  Toeplitz  matrix,  we  determine  the  minimum  necessary  number  of 
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Table  1. 


.13280ill«00 

.22840909<f00 

•59000987-01 

• 105I0259<*>00 

• 102476944>00 
.20997367«00 
.15796515^00 
.49766531^00 
.21688403+00 
.79068109-01 
.36745383400 
.38690715400 
.20062616400 
.79061432400 
.48378490400 
.44248021-01 
.57159514400 
.64039297400 
.23281741400 
.94323532400 
.66032016400 
.62646132-01 
.66972719400 
.80709314400 
.19501743400 
.86397782400 
.64975729400 
.13480165400 
.59166194400 
.61631014400 
.11786686400 
.59562206400 
.45695585400 
.15024535400 
.39399029400 
.65310303400 
.62190854-01 
.30626509400 
.19662871400 
.12695771400 
.20990735400 
.40971760400 
.72522793-01 
.16994681400 
.56496797-01 
.11100784400 
.14634013400 
.23342173400 


*In  row  order. 


Normalized  Optimum  Anq)litude  Coefficients* 


.82793002-01 
.79946042-01 
.135004434>00 
.11804765400 
.12952937400 
.69879676-01 
.20769910400 
.36750016-fOO 
.16740339400 
. 16953667400 
.25321673400 
.24632016400 
.38681863400 
.75863105400 
.23120761400 
.30107630400 
.52915861400 
.51689137400 
.49893874400 
.10000000401 
.29791722400 
.34641811400 
.76620273400 
.70139802400 
.47937338400 
.95812562400 
.33425427400 
.29277998400 
.80448604400 
.72969721400 
.34527393400 
.66302999400 
.30569579400 
.16197544400 
.63671771400 
.50365177400 
.16294750400 
.28356918400 
.21663736400 
.34734902-01 
.37133591400 
.34537617400 
.79662635-01 
.19371152-01 
.13473249400 
.97182458-01 
.15297552400 
.17091950400 


.69682619-01 

.12954030400 

.11604084400 

.13499*^58400 

.79948200-01 

.82793120-01 

.24633li0400 

.25325102400 

.18954971400 

.16736566400 

.36755613400 

.20769469400 

.51691626400 

.52921947400 

.30113671400 

.23113395400 

.75657050400 

.36680651400 

.70143386400 

.76629516400 

.34651732400 

.29782193400 

.99991325400 

.49692030400 

.72973312400 

.80458716400 

.29288662400 

.33415043400 

.95803249400 

.47935589400 

.58367705400 

.63679373400 

.16204624400 

.30562235400 

.66296146400 

.34526428400 

.34538620400 

.37137629400 

.34719955-01 

.21659505400 

.28355505400 

.18294634400 

.17091981400 

.15298356400 

.97175206-01 

.13472130400 

.19381611.01 

.79663160-01 


.20937963400 

.16248319400 

.10518964400 

,53011604-01 

.22640452400 

,13280116400 

.36892712400 

,36746628400 

.73089206-01 

.21686434400 

,49764040400 

,15796425400 

.64043237400 

.57166698400 

.44161563-01 

.46372232400 

.79056159400 

,20082197400 

,80714849400 

.66983428400 

.62609016-01 

.66022555400 

.94316290400 

.23281107400 

,81636749400 

,59177101400 

.13482484400 

.64964522400 

.86390107400 

.19501313400 

.65314620400 

,394n7163400 

.15030165400 

.45887188400 

.59556410400 

,11788674400 

,40973934400 

,20994914400 

.12699706400 

,19659228400 

.30623423400 

,62191825-01 

,23342725400 

.14635226400 

.11101475400 

,565''8338-01 

,16993779400 

.72522303-01 
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Table  2.  Normalized  Optimum  Phase  Coefficients*  ! 


tl2a549l9<»01 

.49387090401 

-.47052112400 

.32125690401 

.a0B36i47«00 

.47209648401 

-.46518051400 

,38068202<f01 

.2825l8574'01 

.29809950401 

.20421724401 

.21630821401  ! 

.21630404-fOl 

.20422223401 

.29810061401 

.2625i920'f0l 

.38066703-»01 

-.46519712400 

.47210894401 

,80839306400 

.32125762>01 

-.47053364400 

.49387190401 

.12854825401 

. 3548969440 1 

.546822994^00 

.33924749401 

.63557684-01 

.34180756401 

.54910171400 

.32653424401 

.16945398400 

.38884690401 

-.22880358400 

.39626531401 

.44318333401 

.44313324401 

.39628450401 

-.22878817400 

.38886137401 

.16950333400 

.32653958401 

.54916155400 

,34181024401 

.63572198-01 

.33924984401 

.54684073400 

.35489718401 

-.17386448400 

.31979584401 

.56545347.01 

.31003526401 

.15524596-01 

.33750664401 

.20219657400 

,31234733401 

.34980986400 

.30915414401 

.24379456400 

-,•57012856-01 

-.57514042-01 

.24393511400 

.30916221401 

,34990495400 

.31235075401 

.20225623400 

.33751057401 

,15543014-01 

.31003644401 

.56550222-01 

.31979665401 

-.17386428400 

.24623131401 

-.34038705400 

.29765801401 

-.13785109400 

.29917011401 

.00000000 

.32453916401 

-.14612001400 

.32253696401 

.10509497400 

.30279628401 

.16217330401 

.16229149401 

.30280263401 

.10523209400 

,32254537401 

-.14610034400 

.32454392401 

.24110079-04 

,29917095401 

-.13784358400 

.29765835401 

-.34039295400 

.24623001401 

.50968453401 

.24662596401 

-.38773412400 

. 28891 855-fOl 

-.28782260400 

.29277679401 

-.25810927-01 

,29196191401 

-.14921746400 

.33210039401 

-.42176142400 

.41287593401 

.41294698401 

-.42180148400 

.33211495401 

-.14919665400 

.29196350401 

-.25782436-01 

.29277723401 

-.28782484400 

.28891654401 

-.38774121400 

.24662427401 

.50968158401 

.13089936401 

-.10508331401 

.24907595401 

-.38131097400 

.27692373401 

-.42415571400 

.29654022401 

-.20916614400 

.27330940401 

.15716442400 

.23393358401 

,47723716400 

.47759348400 

.23390914401 

.15728939400 

,27330796401 

-.20915961400 

.29654120401 

-.42417505400 

,27692279401 

-.38131922400 

.24907413401 

-.10508618401 

,13089531401 

.33542584401 

.14905351401 

-.10180405401 

.26354111401 

-.26036361400 

.24877763401 

-.42787129400 

,31132537401 

-.84244648400 

.31406809401 

.40480070401 

.29940732-t>0l 

.29941207401 

.40468896401 

.31407327401 

-.84256625400 

.31133039401 

-.42788979400 

.24077248401 

-.26036507400 

.26353973401 

-.10180700401 

.14905012401 

.33542486401 

.48997399401 

.31972951401 

.13941491401 

-.85492179400 

.32418062401 

.24888511401 

.20045602401 

.221696054-00 

.14879328400 

-.55421704400 

.51622863401 

.51416783401 

.51416118401 

.51623651401 

-.55423763400 

.14870131400 

.22173604400 

.20045041401 

.24890826401 

,32418243401 

-.65493872400 

.13941286401 

.31972921401 

,48997233401 

*In  row  order. 
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elements  and  set  op  the  complete  Toeplitz  matrix  by  a simple  program 
loop,  then  perform  the  inversion  by  Gauss-Jordan  reduction  subroutine. 

In  reality,  the  individual  directionality  of  the  elements  should 
be  taken  into  consideration.  Incorporating  R|(6,^)  into  the  integrand 
not  only  makes  it  difficult  to  obtain  a closed-form  solution  but  also 
increases  the  computation  time  tremendously.  To  obtain  each  bjj^  through 
numerical  integration  with  element  directionality  will  take  at ''least 
three  times  longer  than  the  omnidirectional  case,  thus  requiring  about  18 
hours  computer  time  to  reach  a solution.  This  time  requirement  almost 
makes  the  solution  impossible.  At  present,  we  use  the  optimum  coeffi- 
cients, obtained  from  the  case  of  omnidirectional  elements,  to  process 
for  an  improvement.  Numerical  results  show  an  inprovement  of  0.28  dB. 
However,  this  answer  remains  to  be  justified  theoretically.  A set  of 
nonoptimum  and  improved  beam  patterns  of  directional  elements  is  shown 
in  figure  4 with  the  associated  DI  values. 


CONCLUSIONS 


Even  though  certain  simplifications  have  been  made  on  the  actual 
physical  conditions  in  order  to  develop  the  SRR  model,  the  model  works 
for  both  passive  emd  active  sonar  systems  in  general.  This  model  deter- 
mines a set  of  coefficients  that  not  only  can  be  used  for  the  maximiza- 
tion of  SRR,  which  leads  to  the  improvement  of  sonar  performance  of 
existing  sonars,  but  also  can  be  used  to  guide  the  design  of  future 
sonars . 

When  modeling  a solution  for  any  conplicated  problem,  it  is  desira- 
ble to  .begin  with  reasonable  simplifications  of  the  actual  problem  to 
make  it  solvable  by  a simple  mathematical  model  with  potentiality  for 
generalization.  When  a mathematical  expression  does  not  have  closed- 
form  solutions,  it  should  be  solvable  by  some  numerical  technique.  It 
is  most  desirable  that  the  solution  be  either  unique  or  optimum.  If  SRR 
can  be  maximized,  the  SRR  model  will  produce  optimum  results. 

It  may  seem  impossible  to  construct  a general  model  to  handle  any 
situation  pertinent  to  a complicated  problem.  Then  it  is  possible  to 
specialize  the  model  for  particular  applications  so  that  advantage  can 
be  taken  of  special  properties.  The  SRR  model  has  been  specialized 
for  a number  of  applications.  Numerical  results  and  speed  show  great 
efficiency  in  this  model. 

In  designing  a computer  model,  it  is  desirable  to  automate  the 
interactions  internally  so  that  the  users  are  not  bothered.  Some  price 
must  be  paid  to  achieve  this  automation.  Whether  or  not  it  is  worth 
the  overhead  to  achieve  the  general  product  has  to  be  weighed  carefully 
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by  the  designer.  The  SRR  model,  as  illustrated  in  the  flowchart  (figure  ] 
2),  is  performed  in  three  stages  with  little  user  intervention  between  | 
stages. 

Modeling  a computerized  mathematical  solution  for  any  problem  is  an  ] 
art.  Besides  accuracy,  speed,  and  mathematical  validity,  a model  must  j 
also  try  to  be  short  and  simple,  as  general  as  possible,  flexible  and  I 
reliable,  and  easy  to  use  in  order  to  be  considered  efficient.  | 
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